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ABSTRACT 

We study non-thermal emissions from cool cores in galaxy clusters. We 
adopted a recent model, in which cosmic-rays (CRs) prevail in the cores and 
stably heat them through CR streaming. The non-thermal emissions come from 
the interaction between CR protons and intracluster medium (ICM). Comparison 
between the theoretical predictions and radio observations shows that the overall 
CR spectra must be steep, and most of the CRs in the cores are low-energy CRs. 
Assuming that the CRs are injected through AGN activities, we study the nature 
of the shocks that are responsible for the CR acceleration. The steep CR spectra 
are likely to reflect the fact that the shocks travel in hot ICM with fairly small 
Much numbers. We also study the dependence on the CR streaming velocity. 
The results indicate that synchrotron emissions from secondary electrons should 
be observed as radio mini-halos in the cores. In particular, low-frequency obser- 
vations (e.g. LOFAR) are promising. On the other hand, the steepness of the 
spectra makes it difficult to detect non-thermal X-ray and gamma-ray emissions 
from the cores. The low-energy CRs may be heating optical filaments observed 
in the cores. 

Subject headings: cosmic rays — galaxies: clusters: general — galaxies: clusters: 
intracluster medium — radiation mechanisms: nonthermal 
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Introduction 



Clusters of galaxies are filled with hot X-ray gas or intracluster medium (ICM) with tem- 
peratures of ~ 2-10 keV. While the radiative cooling time of the ICM is longer than the age of 
the Universe in most of the region in a cluster, the exception is the core, which is r < 100 kpc 
from the cluster center. If there is no heating source, the ICM in the core cools and a flow 
toward the cluster center should develop (a cooling flow). However, X-ray observations 



should be heated by some unknown sources (e.g. 


Ikebe et al. 


1997 




Makishima et al. 


2001; 


Peterson et al. 


2001; 


Tamura et al. 


2001; 


Kaastra et al. 


2001 


: : 


Matsushita et al. 


2002 


). Since 



active galactic nucle i (AGNs) are often found in the cores, they are often thought to be the 



heating sources (e.g. IChurazov et al.ll200ll ; iQuilis. Bower. &; Baloghll200ll ; iBriiggen fc Kaiser 



tion between AGNs and the ambient ICM (e.g. 


Fabian et al. 


2000; 


McNamara et al. 


2000 




Blanton et al.lboOl: 


McNamara et al.lboOl 


Mazzotta et al. 


200 


2;lFuiita et al. 


2002 


Johnstone et al. 


2002; 


Kempner. Sarazin. & Ricker 


20021 


Takizawa et al. 


2003; 1 


Jujita et al. 


2004b|). 


However, 



even if AGNs can produce enough energy to heat the core, the energy must deliberately be 
transported to the surrounding ICM in the core. For example, conventional mechanical heat- 
i ng such as the dissipat i on of weak shocks and sound waves often c ause thermal instabilities 
( jFujita fc Suzuki! 120051 ; iMathews. Faltenbacher. fc Brighentil 120061 ). Therefore, turbulence 
may essentially be required to hold the instabilities for such heating mechanisms. 



Cosmic-rays (CRs) may be anot 



Tucker &; Rosner 



2004 



1983 



Pfrommer et al. 



Rephaeli 



2007 



198 



7: 



rer channel of transporting energy to the ICM (e.g. 



Jubelgas et al. 



as an < 



Rephaeli &; Silklll995l ; IColafrancesco. Par, fc De Rujula 



2008). Especially, CR streaming has been studied 



energy transport mechanism ( Rephaeli 19 79 ; Bohringer fe Morfill 1988 ; Loewenstein. Zweibel. &; Beg' 



199 It iGuo fc Ohll2008h . In this mechanism, CRs streaming in the ICM excites Alfven waves. 
The CRs interact and move outwards with the waves. The PdV work done by the CRs 
effectively heat the ICM. Recently, using numeric al simulations, we sh owed that the CR 
streaming can stably heat the core for a long time ( jFujita fc Ohirall201ll hereafter Paper I). 
The reason of the stability is that the CR pressure is insensitive to changes in the ICM 
and that the density dependence of the heating term is similar to that of radiative cooling. 
Moreover, CRs can prevail in the entire core and the heating is not localized around the 
source. The C Rs may be provided in the core not only by AGNs but also through pumping 
by turbulence ( jEnBlin et al.l 120111 ). 



In this paper, we study the non-thermal emission from the CRs that heat cool cores and 
the AGN activities that are responsible for the acceleration of the CRs. It is to be noted 
that non-thermal emissions from CR protons accelerated by AGNs in the cores have been 



-3 - 



studied by iFujita et al.l (120071 ). However, they studied CR acceleration associated with a 
single AGN burst with an extremely large energy, and they did not consider the heating of 
the ICM by CR streaming. This paper is organized as follows. In § 121 we explain our models 
on core heating and AGN activities that are responsible for the generation of CRs. In § [21 
we present the results of our calculations and compare them with observations. In § HI we 
discuss the implications of our results, and §[5] is devoted to conclusions. We refer to protons 
as CRs unless otherwise mentioned. 



2. Models 



2.1. Cosmic-Ray Distributions 



In Paper I, we studied heating of a cool core by CRs injected through the activities of 
the central AGN. The CRs travel with Alfven waves in the ICM. They amplify the waves, 
which heat the surrounding ICM. In this subsection, we briefly summarize the models to 
obtain CR and ICM distributions. 



For simplicity, we assumed that the cluster is spherically symmetric. The flow equations 



are 



o 



dt r 2 dr r 2 dr 



Pg + P c + P B ) 



(1) 

(2) 



de„ 1 d . 2 , 



1 d 



-P^jr^u) + — — 



' 2 m ar' 



-n 2 e A{T) + H st + H coll 



— + — — (r 2 ue ) 
dt r 2 dr c 



r 2 D(p) 



de c 
dr 



Tloss S c 



(3) 
(4) 



where p is the gas density, u is the gas velocity, P g is the gas pressure, P c is the CR pressure, 
Pb is the magnetic pressure, G is the gravitational constant, M(r) is the gravitational mass 
within the radius r, k(T) is the coefficient for thermal conduction and T is the temperature, 
n e is the electron density, A is the cooling function, H st is the heating by CR streaming, H co u 
is the heating by Coulomb and hadronic collisions, u is the CR transport velocity, D(p) is 
the diffusion coefficient for CRs averaged over the CR spectrum, T\ oss is the energy loss by 
Coulomb and hadronic collisions, and S c is the source term of CRs. Energy densities of the 
gas and the CRs are respectively defined as e g = Pg/i'jg — 1) and e c = P c /{lc — 1), where 
7 S = 5/3 and 7 C = 4/3. In this paper, we do not treat models with thermal conduction, and 
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thus k — 0. The terms for radiative cooling A, Coulomb collisions H co u, hadronic collisions 
if co ii, diffusion D(P), and the energy loss r loss are the same as those in Paper I. The source 
term of CRs is given by S c oc Lagn, where Lagn is the energy injection rate from the AGN. 
We assume that Lagn = eMc 2 , where e is the parameter, M is the inflow rate of the gas 
toward the AGN, and c is the speed of light. 



The CR transport velocity in equation (j3J) is given by u = u + va, where va = B/^Airp 
is the Alfven velocity for a magnetic field B, which evolves as B oc p 2 / 3 . The initial magnetic 
field at the cluster center is B = 10 fxG. The wave energy U A = 5B 2 /(An), where SB is the 
magnetic field fluctuation, is amplified by the PdV work done by the CRs on Alfven waves: 



dU A | 1 
dt r 2 



d_ 

dr 



r 2 U A [ -u + v A 



u 



8_Ua 
dr 2 



- v A 



dPc 
dr 



-H, 



St 



(5) 



(IMcKenzie &: Volk!ll982l ; iBohringer fc Morfilll l 19881 ) . This equation is more correct than tha ; 
we adopted in Paper I (equation [6] in that paper), because it is based on wave energy conser- 
vation. However, the results are not affected by this change of the equation (see §|3j). After 
the wave energy increases to Ua ~ Um, where Um is the energy of the backgrou nd magnetic 
field, the waves are expe cted to heat ICM through non-linear damping (e.g. lOhira et al. 
20091 ; iGargate et al.l 120101 ) . Thus, we give the heating term for CR streaming by 



dP r 



dr 



flVolk et al.lll984J : iKang fc Jonesll2006f ). We simply give T 
after Ua reaches Um- 



(6) 



Ua/Um for Ua < Um and T = 1 



2.2. Non-Thermal Emissions 

Although by solving equations presented in § 12.11 we can obtain the profile of the ICM 
and that of the CR pressure P c {r) required to heat the core effectively (Paper I), we do not 
have information on the energy spectrum of the CRs. Thus, we need to specify the spectrum 
of the CRs to calculate the non-thermal emissions from the CRs. 

We assume that the central AGN drives outgoing shock waves and form cocoons or 
bubbles inside them. In the following, we show a description of their evolution (position, 
velocity, Mach number as functions of time). The shocks should inject CRs with varying 
efficiencies and spectra. At some moment this injection will be maximal (actually, this 
moment differs for different CR energy ranges). We only consider the efficiency and Mach 
number at this moment and fix these numbers by requesting them to reproduce the observed 
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radio emission for the sake of simplicity, although there would be a more physical approach 
to calculate the injection evolution and the full injected CR spectrum, assuming the AGN 
energy release, timescale, and initial cocoon radius. 

The CRs are accumulated in the core through the AGN activities. In Paper I, we 
studied continuous CR injection as a time average, although the supply of the CRs may 
be intermittent. Each activity of the AGN is approximated by an instantaneous explosion. 
Thus, the shock expands in the ICM like a supernova remnant in the Galaxy, and the shock 
velocity depends on the energy input from the AGN. 

In Paper I, we obtained the profiles of the ICM density p(r), the temperature T(r), and 
the magnetic field B(r) at a given time. We approximate the density profile of the ICM by 
a power-law: 

PicmW = Pm(r/r in )~ u , (7) 

wher e p\ n is the ICM density a t the inner boundary (r- m = 5 kpc). Using a shell approximation 
(e.g. 



Ostriker fc McKedll988f ). the radius of the shock can be written as 



where 



R s 



id 



E n 



Pin r in 



1/(5-0)) 



+2/(5-0,) 



i_(T 9 + l) 2 (7 s -l)(3-o;) 

47T 



l/(5-w) 



(9) 



9 79 -3-u;(7 9 + l) 

E a is the energy released by the AGN, and t a is the time elapsed since the last energy input 
from the AGN. The velocity of the shock is given by 



V. 



dR s 
dt a 



(10) 



The Mach number of the shock is given by M s = V s /c s (R s ), where c s is the sound velocity. 
Since we know the profile of the ICM temperature T(r), we can construct the profile of the 
sound velocity c s (r). Therefore, if M s and E a are given, the shock radius R s , velocity V s , 
and the time t a that satisfy equations (JH]) and ( flOl can be specified. 

In reality, the spectrum of accelerated CRs at the shock may change during the expansion 
of the cocoon. Probably, the spectrum is flat, when the cocoon is young, and the shock 
velocity and the Mach number are large. Then it gradually steepens as the Mach number 
decreases, and CR acceleration ceases when the Mach number approaches M s ~ 1. However, 
we consider a typical Mach number M st around which most CRs are accelerated. In other 
words, we consider a typical spectrum of CRs that are accelerated when the injection of CRs 
becomes maximal. We treat M st and E a as parameters. The shock radius, velocity, and age 



-6 - 



when M s = M st are R s = R st , V s = V st , and t a = t at , respectively. We also assume that the 
spectrum of CRs that are just accelerated at r ~ R st has a form of 

N(p, r) oc p- x e~ p/pxam ' , (11) 

where p is the CR momentum, x is the index, and p max is the cutoff momentum of the CRs. 
Since we already know CR pressure P c {r), the normalization of relation (II ip is determined 
by the relation 

c f°° p 2 N(p,r) 

r-cKH = 

where m is the proton mass. 



P c (r) = - \ - K = ' dp , (12) 

3 Jpndn VP 2 + m2 ° 2 



T he index is given by x = ( r& + 2)/(r& — 1), where r& is the compression ratio of the 



shock (IBlandford fc Eichlerl 1 1 98 7t ) . which is given by 



n - {l9 + l)Ml (13) 

Since the cooling of CR protons is not effective, the maximum energy of protons correspond- 
ing to p max is determined by the age of the shock and is represented by 



where is the downstr eam magnetic field at r = R st and is given by = r^B (jYamazaki et al 



20061 ; iFujita et al.ll2007l ). That is, the background magnetic field B is amplified by the com- 
pression ratio Tb (equation [T5]). Although some particles may be accelerated to higher 
energies when the expansion velocity of the cocoon was larger, their contribution to the 
overall spectrum is expected to be small. 

The CRs injected at r ~ R st propagate in the ICM with Alfven waves. Although 
adiabatic cooling may change E max , it does not change the index x in relation (iTTj) . Moreover, 
the results in §[3]show that the CR spectra must be steep. Thus, the results are not sensitive 
to the value of E max . Therefore, we do not consider the adiabatic cooling for E ma:K and adopt 
the relations f llip and ffT2"j) at any radius r, although the adiabatic cooling was considered 
when we calculated P c in Paper I. Since we expect that thermal protons with higher energies 
are accelerated as CRs, we assume that the minimum momentum of the CRs is p m i n = 4mc s d, 
where c sa > is the sound velocity of the ICM at the downstream of the shock at r = R st , which 
is obtained from the Rankine-Hugoniot relations for given c s (R st ) and M st : 



Csd - ° S (7, + l)Mst ' ( } 
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In this way we have the CR spectrum at each radius for given M st and E a . 



For a given CR proton spectrum, we calculate radiation from them. Non-thermal emis- 
sion origina ted from CR protons in the central region of clusters have been studied by several 
groups (e.g. lMiniatill2003l; iPfrommer fc En6rinll2004[ iKeshet fc Loebll2010l ). In this paper, we 



adopt the model of lFuiita et all J2OO9L in which they calculated non-thermal emissions from 



supernova remnants. We consider the synchrotron, bremsstrahlung, and inverse Compton 
(IC) emissions from secondary electrons created through the decay of charged pions that are 
generated through proton-proton collisions. IC emissions are created by electrons that scat- 
ter Cosmic Microwave Background (CMB) photons. We also consider 7r°-decay gamma-rays 
through proton-proton collisions. We do not consider emissions from primary electrons that 
are directly accelerated at the shock, because we did not calculate the distribution of the 
primary electrons in Paper I. Because of the short cooling time of electr ons, emissions from 
primary electrons will disappear soon after their acceleration is finished (IFujita et al.l 120071 ). 



The photon spectra are calculated based on the radiation models of iFang fc Zhang 



(120071 ) . For the production of secondary electrons and 7r°-d ecay gamma-ray photons t hrough 
proton-proton interactions, we use the code provided by iKarlsson fc Kamad (120081 ). The 
spectrum of the secondary electrons are given by N e (E e ) = t coo \ te (E e ) Q e (E e ), where E e is 
the electron energy, t coo i,e is the cooling time of an electron, and Q e is the production rate 
of the secondary electrons. For the cooling, we include synchrotron cooling, IC scattering, 
Bremsstrahlung, and Coulomb loss. 



3. Results 

Since we replaced the equation for the wave energy Ua (equation [6] in Paper I) with 
equation ([5]), we recalculate the distributions of the ICM and CRs and show them in Figures[T] 
and [2j The cluster is initially isothermal with P c = 0. The input parameters are the same 
as those of Model LCR0 in Paper I, and we simply refer to this model as LCR0 again. The 
gravitational potential adopted in this model is that of the Perseus cluster. The efficiency 
of AGN energy input is e = 2.5 x 10 -4 . The results are almost identical to those in Paper I 
(Figures 2 and 4 in that paper). This is because T rapidly approaches one after Lagn 
increases regardless of the equation we adopted. The ICM temperature outside the core 
is ~ 7 keV. The heating by CR streaming and the radiative cooling are well-balanced at 
t > 4 Gyr. 

Figure [3] shows the spectra of a region including the entire core (r < 1 Mpc) at t — 9 Gyr 
for Model LCR0. The slope of the ICM density profile is assumed to be u = 1, which is a good 
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approximation for r < 70 kpc (Figure [Qd). The distance to the cluster is 78.4 Mpc, which is 
the one for the Perseus cluster. In this figure, we take M st = 2.1 and E a — 1 x 10 60 erg s -1 ; 
we first give E a , and then adjust M ?t in order to be consistent with radio obser vations for 
the mini-halo in the Perseus cluster ( 1SijbringHl993t iGitti. Brunetti. fc Settill2002l ). Since the 
Mach umber M st is fairly small, the CR spectrum is steep (x = 3.2). The maximum energy 
of the CRs is E mSuX = 1.5 x 10 5 TeV, the radius and age of the shock are R st = 22 kpc 
and t at = 6.0 x 10 6 yr, respectively. The spectrum of thermal Bremsstrahlung is shown for 
comparison. 

The slope of the synchrotron and IC scattering spectra at the higher energy side can be 
explained by a simple calculation. The slope of the energy spectrum of secondary electrons 
are the same as that of protons (x = 3.2), if radiative cooling is not effective. However, 
cooling by sync hrotron radiat ion and IC scattering increases the slope by one and it becomes 



x 



4.2 (e.g. ISarazinl Il999l ). The spectral indices of the synchrotron emissio n and IC 



scattering are represented by a = (x' — l)/2 
consistent with those in Figure |3] (f v oc v~ a ) 



1.6 (e.g. iRybicki fc Lightmanlll979l ). which is 



We found that the results for M st = 2.1 and E a — 1 x 10 61 erg s _1 are almost the same 
as those for M st = 2.1 and E a = 1 x 10 60 erg s -1 . In the former case, the shock radius and 
age are R s t = 59 kpc and t at = 1.4 x 10 7 yr, respectively. Although the maximum energy of 
the CRs is increased (-E max = 2.7 x 10 5 TeV), the steep CR spectrum or the large x obscures 
the effect. This means that the radiation from the cool core is insensitive to the strength 
of an AGN activity (E a ) for a given P c (r). Figure H] shows the spectra when M s = 1.8 and 
4.0 for E a = 1 x 10 60 erg s _1 . The indices are x = 3.8 and 2.3, respectively. Compared with 



the results of M s = 2.1 and E a = 1 x 10 60 erg s~\ the non-thermal emissions are weaker 
(stronger) when M s is smaller (larger), and the synchrotron radio emission is inconsistent 
with the observations. The results are very sensitive to the value of M st for a given P c (r). 
Basically, changing E a and M st correspond to changing E mgx and x, respectively. 

Figure |5] shows the surface brightness profiles for Model LCRO at t — 9 Gyr with 
M st = 2.1 and E a = 1 x 10 60 ergs _1 . The model generally reproduces the surface brightness 
profile observed in the radio band, although we did not intend to reproduce that when we 
calculated Model LCRO in Paper I. In that figure, the surface brightness rapidly increases 
toward the cluster center for the synchrotron radio emissions because of the increase of the 
magnetic fields toward the cluster center (B(r) oc p(r) 2 / 3 ). On the other hand, the pro- 
file for the IC missions is relatively flat because electrons scatter CMB photons, which are 
uniformly distributed. The size of the region with high surface brightness is regulated by 
radiative cooling, because radiative cooling increases the ICM density and makes a cool 
core. On the other hand, CRs can fill the entire core with fast Alfven waves (Paper I). 
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Proton-proton interactions are effective in such a high-density region. The surface bright- 
ness for thermal Bremsstrahlung slightly decreases at the cluster center, because the ICM 
temperature decreases there (Figure Hk). 

We also study a less massive cluster. Figure |6] shows the spectra of the entire core 
(r < 1 Mpc) at t — 9 Gyr calculated using parameters of Model SCRO in Paper I. For this 
model, we adopted the observed gravitational potential of the Virgo cluster. The efficiency of 
AGN energy input is e = 1 x 10~ 4 . Although we recalculated the ICM and CR distributions, 
they are almost identical to those calculated in Paper I. The ICM temperature outside the 
core is ~ 2 keV. We take M st = 2.1 (x = 3.2) and E a = 1 x 10 59 erg s -1 . The distance 
to the cluster is set to be 16 Mpc. We take u = 0.7, which is an good approximation 
for r < 50 kpc (Figure 11 in Paper I). The shock radius and age are R st = 21 kpc and 
t at = 6.3 x 10 6 yr, respectively. The maximum energy of CRs is E max = 3.3 x 10 4 TeV. 
The gamma -ray flux is much small er than the upper limits for the Virgo cluster obtained 
with Fermi (lAckermann et al.l 120101 ) . The luminosity is sensitive to M st but not to E a . The 
surface brightness profiles for this model are shown in Figure [7J The surface brightness is 
smaller than that in Figure 



4. Discussion 



We have studied the non-thermal spectra of cool cores heated by CR streaming. The 
results indicate that the Mach number of the shock that accelerate CRs must be small (~ 2) 
to be consistent with radio observations at least for the Perseus cluster. We think that this 
is reasonable because the temperature and the sound velocity of the ICM is large and thus 
it is difficult for the cocoon shock to have a large Mach number. The small Mach number 
means that the CR spectrum must be steep. In Paper I, we did not specify the injection 
mechanism of CRs. We emphasize that even if CRs are injected by anything other than the 
cocoon, the spectrum must be steep for the given P c . 



Recently, lEnBlin et al.l ( 1201 ll ) indicated that the CR streaming velocity may be much 
larger than the Alfven velocity va in the hot ICM. This is because in high-/? plasma, where 
(3 is the ratio of thermal to magnetic energy, waves may suffer strong resonant damping by 
thermal protons. In t his case, the sound velocity c, may be appropriate as the streaming 
velocity instead of va (IHolman et al.lll979t lEnBlin et al.l 1201 ll ). Thus, we simply replace va 
by c s in equations (TJJ and ([5]) and see what would happen. Figures [8] and [9] show the profiles 
of the ICM and CRs for the parameters of Model LCR0 except for the larger streaming 
velocity c s ; we refer to this model as Model LCRs. The ICM is stably heated by the CR 
streaming even in this model and the evolution of M is not much different from that in Model 
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LCRO. Compared with Figures [21 the fraction of CR pressure is small in the central region 
because of the larger streaming velocity and the escape of CRs. Since the ICM temperature 
is an increase function of radius (Figure [S]), the sound velocity or the streaming velocity is 
also an increasing function. Thus, P c /P g tends to decrease outward fairly rapidly. Figures [TU] 
and [TT1 are the same as Figures |3] and |5] but for Model LCRs. If we assume M st = 2.1 as is 
the case of Model LCRO, non-thermal luminosities in Model LCRs are smaller than those 
in Model LCRO, because more CRs have escaped from the core with the high ICM density. 
Thus, we increase the Mach number and set it to be M st = 2.4. We present the spectra 
and surface brightness in Figures [10] and [TTJ The synchrotron spectrum and the surface 
brightness are consistent with the observations. 

Regardless of the streaming velocity, the CR spectra in the cores must be steep, because 
if not, the luminosities are too large (Figure Sb); this is in consistent with the small number 
of clusters in which radio mini-halos have been observed (IGovoni et al.ll2009l ) and the non- 
detection of gamma-rays from clusters. Because of the steep spectra, future observations 
in the low-frequency radio band would be useful. Thus, cool cores would be promising 
targets for radio telescopes such as LOFAR. The number of mini-halos may increase as the 
sensitivities of radio telescopes are improved. In our model, we assumed that CRs are mostly 
accelerated when the Mach number of the shock is M s ~ M st . For real clusters, however, we 
expect that the Mach number M s decreases and that the CR spectrum at the shock steepens 
during the expansion of a cocoon. Thus, we expect that the spectral index should increase 



outwards in the cluster, which has actually been observed in the radio band (ISijbring||l993 



Gitti et al.ll2002l ). although the interferometric nature of these measurements might result in 
smaller radio halos at higher frequencies (missing zero spacing problem). On the other hand, 
observations in other bands would be difficult in the near future (Figures [3] and |6]). In the 
X-ray band, IC emissions should be observed (Figures [3] and [6]). However, thermal emissions 
from cool cores are very bright, wh ich makes it difficult for the non- thermal emission to be 
detected. For the Perseus cluster, ISanders. Fabian, fc Dunnl (120051 ) claimed the detection 



of non-thermal emission with a flux of 6.3 x 10 11 erg cm 2 s 1 between 2 and 10 keV 



However, the detection was not confirmed by later observations (IMolendi fc Gastaldelldl2009 



Eckert fc Paltanill2009l ). Even with hard X-ray telescopes that will be launched in the near 
future such as NuSTAR and Astro-H, the detection may be difficult because of the low 
surface brightness (Figures [5] and [7]). For the detection in the gamma-ray band, good angular 
resolutions as well as sensit ivities are required, because gamma-r ays could also be emitted 
from the central AGNs (e.g. lAbdo et al.ll2009l ; iKataoka et al.ll2010l ). which must be resolved. 



The steep CR spectra mean that most of the CRs in cool cores have low energies. Thus, 
indirect studies may be usef ul. For example, op t ical filaments obs erved in cool cores may 
be heated by those CRs fsee iFerland et al.lbood : bavet et al.lhoioh . We note that the CR 
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heating is locally unstable, and that the filaments could be created through local thermal 
instabilities (Paper I). Moreover, our model does not require turbulence for stable heating. 
Thus, cool cores in which turbulence is not developing may be observed with detectors having 
high spectral resolutions such as Astro-H, while the detection of turbulence does not deny 
our model. Although we did not consider prim ary electrons, they ma y be accelerated at 
shocks in cores in spite of the low Mach numbers (jMatsukiyo et al.ll201ll ). and the emissions 
from them may be observed in some clusters. 

Finally, we caution the reader that we did not consider energy- dependent diffusion 
of CRs, because we do not k now the actual di ffusion coefficient in the ICM, especially 
away from the shock front (see lFujita et al.ll201ll ). If CRs with higher energies escape from 
the core faster than those with lower energies, the energy spectrum could be steep (e.g. 



Fujita et al.ll2009t lOhira. Murase. fc Yamazakill2010l 1201 if ). Moreover, we did not include 



the contribution o f gamma-rays from CRs accelerated at cosmological shocks and those from 
dark- matters (e.g. iTotanil 120041 : iPinzke fc Pfrommerll2010l ; iPinzke. Pfrommer. fc Bergstrom 

201 ih . 



5. Conclusions 



We have investigated non-thermal emissions from cool cores heated by CRs. For the 
distributions of CRs, we used the model in which the cores are stably heated by CR stream- 
ing. CR protons interact with ICM protons and produce secondary electrons and 7r°-decay 
gamma-rays. We found that the CR spectra must be steep in order to be consistent with 
observations of a radio mini-halo. The steep spectra reflect the fact that the CRs are ac- 
celerated at shocks with low Mach numbers (~ 2) in hot ICM. We have also studied the 
dependence on the CR streaming velocity and found that the stronger shocks are required 
to be consistent with the observations for the larger CR streaming velocity. Since most of 
the CRs in cores have low energies, synchrotron emissions from them should be observed in 
low-frequency radio bands. Thus, the number of clusters that have radio mini-halos would 
increase as the sensitivities of radio telescopes increase. On the other hand, the detection in 
other bands such as the X-ray and gamma-ray bands would be difficult in the near future. 
The low-energy CRs could be studied by observing optical filaments that are often found in 
cool cores. 



We thank the referee for useful comments. This work was supported by KAKENHI 
(Y. F.: 23540308, Y. O.: 21684014). 
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Fig. 2.— Profiles of the ratios P c /P g (solid) and P B /P g (dotted) at t = 9 Gyr for 
Model LCRO. 
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thermal bremsstrahlung (dashed line) are of the secondary electrons. The 7r° decay gamma- 
rays are shown by the two-dot-dashed line. The thermal bremsstrahlung is shown by the 
dot-dashed line. Observations are for the Perseus cluster. Radio observations are shown 
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Fig. 4.— Same as Fig. [3] but for (a) M s = 1.8 and (b) 4.0. 
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Fig. 5. — Surface brightness of non-thermal and thermal emissions calculated based on 
Model LCR0 with M st = 2.1 and E a = 1 x lO^ergs" 1 . The synchrotron radiation (327 MHz; 
dotted line), IC scattering (20 keV; solid line), 7r° decay gamma-rays (1 GeV; two-dot-dashed 
line), and thermal Bremsstrahlung (20 keV; dot-dashed line) are shown. Radio observation s 
for the mini-halo in the Perseus cluster are shown by dots (jGitti. Brunetti. &: Settill2003l ). 
The vertical errors include the deviations from the spherical symmetry. 
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Fig. 6. - — Same as Fig. [3] but for Model SCRO. Gamma-ray upper limits are shown by an 
arrow (lAckermann et al.ll2010l ). 
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Fig. 7.— Same as Fig. M but for Model SCRO. 
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Fig. 8- 



Same as Fig. [T]but for the Model LCRs. 
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- Same as Fig. EJ but for Model LCRs and M st = 2.4. 
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Fig. 11- 



Same as Fig. but for Model LCRs and M st = 2.4. 



